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Abstract 

A consistent microscopic diagrammatic approach is applied for the first time to 
the calculation of the nucleon emission spectra in the non-mesonic weak decay of 
A-hyper nuclei. We adopt a nuclear matter formalism extended to finite nuclei via 
the local density approximation, a one-meson exchange weak transition potential 
and a Bonn nucleon-nucleon strong potential. Ground state correlations and final 
state interactions, at second order in the nucleon-nucleon interaction, are introduced 
on the same footing for all the isospin channels of one- and two-nucleon induced 
decays. Single and double-coincidence nucleon spectra are predicted for j^C and 
compared with recent KEK and FINUDA data. The key role played by quantum 
interference terms allows us to improve the predictions obtained with intranuclear 
cascade codes. Discrepancies with data remain for proton emission. 
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Being the only source of information available on strangeness-changing baryon 
interactions, the non-mesonic weak decay of A-hypernuclei has attracted con- 
siderable interest and experienced great advances in the last years [Tf2f3] . 
These hadronic weak interactions, whose determination requires the solution 
of complex many-body problems also involving strong interaction physics, 
are important inputs, for instance, when investigating the thermal evolu- 
tion and the stability of compact stars [3]. The non-mesonic decay width, 
Tnm = Ti + r 2 , is built up from one- (IN) and two-nucleon induced (2N) 
decays, T 1 = T n + T p and T 2 = T nn + Y np + T pp , with the isospin compo- 
nents given by T N = r(AiV ->• nN) and T NN , = T(ANN' ->• nNN'), with 
N, N' = n or p. 
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As noted in [5] , from a quantum-mechanical point of view the only observable 
in non-mesonic weak decay are the rate T^m and the spectra of the emitted 
nucleons. Each one of the elementary non-mesonic decays occurs in the nu- 
clear environment, thus final state interactions (FSI) can modify the quantum 
numbers of the nucleons produced in such weak processes and new, secondary 
nucleons are emitted as well: this prevents the measurement of any of the 
non-mesonic partial decay rates. 

Although the long-standing puzzles on T n /T p and the asymmetry parameter 
have been solved recently [2], a detailed knowledge of the non-mesonic decay 
mechanisms is still missing, as demonstrated by the persistent discrepancies 
between theory and experiment on the nucleon spectra [6]. 

In this Letter we present for the first time a consistent microscopic calcula- 
tion of the nucleon emission spectra. A nuclear matter formalism is adopted 
and results for single and double-coincidence nucleon spectra are reported 
for within the local density approximation. Previous investigations [7] 
demonstrated the negligible effect of ring and RPA microscopic contributions 
on nucleon spectra. Ground state correlations (GSC) and FSI contributions 
are introduced here at second order in the nucleon-nucleon interaction for the 
whole set of IN and 2N isospin decay channels. The weak transition potential 
yAN^NN con t a i ns the exchange of the full set of mesons of the pseudoscalar (it, 
r], K) and vector octets (p, u, K*), with strong coupling constants and cut-off 
parameters deduced from the Nijmegen soft-core interaction NSC97f [S]; for 
the nucleon-nucleon interaction y NN ^ NN we adopt the Bonn potential (with 
the exchange of it, p, a and u mesons) [9]. Being fully quantum-mechanical, 
the present approach is expected to produce more reliable results for FSI than 
those based on the (semi-classical) nucleon rescattering given by intranuclear 
cascade (INC) models [TUlfTT| . 

The many-body terms we consider are represented by the in-medium A self- 
energy Feynman diagrams of Fig. [TJ where we limit to two- and three-nucleon 
emission. Diagrams D and E contribute to IN decays. By considering all the 
possible time-order ings of the Feynman diagrams at second order in y NN ^ NN , 
one obtains Goldstone diagram contributions to 2N and FSI-induced decays 
as well as quantum interference terms (QIT). FSI Goldstone diagrams have at 
least one y NN ^ NN acting after y AN ^ NN . QIT always involve a FSI and are 
of two kinds: FSI amplitudes interfere with both IN (1A^-FSI QIT) and 2N 
decay amplitudes (2A-FSI QIT). We shall use the term 'plain FSI' to indicate 
a FSI which is not a QIT. For example, from the pp Feynman diagram of Fig. [T] 
one obtains the Goldstone diagrams of Fig. [2J diagram (a) provides a plain FSI 
contribution when a cut on 3p2h states is considered (this cut is associated to 
three-nucleon emission), while for cuts on 2plh states (two-nucleon emission) 
one has 1A-FSI QIT; (b) is a 2N decay contribution, introduced by GSC; (ci) 
and (c 2 ) are QIT with two possible final states, 2plh, contributing to 1A-FSI 
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Fig. 1. The set of Feynman diagrams considered in this work for the in-medium A 
self-energy. D and E are the direct and exchange terms of the IN decay channel, 
the remaining ones contribute to 2N and FSI-induced decays. The dashed and 
wavy lines stand for the weak and strong potentials, y& N -* NN and y NN ~* NN ; 
respectively. 

QIT, and 3p2h, leaving to 2iV-FSI QIT (incidentally, they are vanishing since 
Pi = hi and then one of the two requirements on the Fermi level, pi > kp and 
hi < kp, cannot be fulfilled). 

(«) (b) (ci) (c 2 ) 

L I _ L_ I. 





Fig. 2. The Goldstone diagrams corresponding to the pp Feynman diagram of Fig. [TJ 
Note that the QIT (ci) and (02) are vanishing since pi = hi and then one of the two 
requirements on the Fermi level, Pi > kp and hi < kp, cannot be fulfilled. See text 
for details. 



From this analysis it is clear that a single Goldstone diagram can contribute to 
a plain FSI and to a QIT, depending on the final physical state one considers: 
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this point is important to have an unambiguous interpretation of the QIT. A 
cut in a Goldstone diagram that leads to two symmetric pieces is interpreted as 
the square of a transition amplitude: a non-QIT which contains the interaction 
yNN-tNN can re p resen t either a 2N decay or a plain FSI-induced decay. At 
variance, when the cut divides the diagram into two different amplitudes, one 
has a QIT. 2N and plain FSI contributions are positive-definite, while QIT 
can be either positive or negative. This trivial fact has important consequences 
that will be discussed later. 

Regardless of the set of many-body contributions one considers, the Goldstone 
diagram technique allows us to write the total number of nucleons and nucleon 
pairs emitted in the non-mesonic decay as follows [7j: 

N n = 2T n + T p + 3T nn + 2T np + T pp + Nf ( n )^i,f > 

i,f 

Np = T p + T np + 2f p P + ^Nf ( p ) Tij , 

ij 

■N n n -Tri -\- 3r nrl -\- -|- ^ K Nf (jm) I" \ t f , 

hf 

N n p Tp ~\~ 2r n p -)- 2Tpp -|- ^ Nj (np) ^i,f j 

Np P = f pp + N f (pp) f ij , ( 1 ) 

i,f 

where a normalization per non-mesonic decay is used (f = r/r NM ). Single 
and double coincidence nucleon spectra are obtained by constraining the eval- 
uation of each f to certain intervals in energy, opening angle, etc. The Tjv's 
(rViv's) are the IN (2N) decay rates, while the remaining terms containing 
the functions Tij represent FSI Goldstone diagrams (for instance, diagrams 
(a), (ci) and (c 2 ) of Fig. [2]). The index i in Tij is used to label the FSI Gold- 
stone diagrams obtained from the Feynman diagrams of Fig. [T] at second order 
in v NN ~^ NN , while / denotes the final physical states of the Goldstone dia- 
gram and in the present case can take the values / = nN (cut on 2plh states) 
and nNN' (cut on 3p2h states). Finally, ^V/(jv) (^f(NN')) is the number of 
nucleons of the type iV (of iViV' pairs) contained in the multinucleon state /. 

Note that Eqs. ([1]) contain five observables, N^- and iVVv , and five unobserv- 
able decay widths, T N and T NN /. The unobservable character of the decay 
widths is due to the presence of FSI, which we know to be important. There- 
fore, one cannot invert the above relations to obtain the decay widths in terms 
of the observed spectra. The experimental values of the partial widths can be 
obtained from the measured spectra only if a model is used for deconvoluting 
the FSI effects. 

To illustrate an important property of FSI, let us consider the diagram of 
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Fig. (we assign to it the index i = 1) with / = nNN 1 (cut on 3p2h 
states): it is represented by the function Ti^nn 1 , which is a plain FSI-induced 
contribution. The diagram Ufa) can be cut also on 2plh states: in this case 
it contributes to the liV-FSI QIT F li7l jy From sum rule considerations on 
Goldstone diagrams [12] it follows that 

^ f l, n JV + X! ^l,nNN' = 0. (2) 
N=n,p N,N'=n,p 

In general, each FSI A self-energy diagram can in no way contribute to the 
total non-mesonic rate: 

2J ^i,nN + J! ^i,nNN' + ^i,nNN'N" + • • • = 0, (3) 

N=n,p N,N'=n,p N,N',N"=n,p 

for any value of i, while in principle it can affect the IN, 2N, etc, decay rates 
separately (see however the discussion of [5] on the opportunity not to include 
these FSI terms in the definition of the unobservable partial rates). Note that 
the above sum rule is valid for the nucleon spectra integrated over the whole 
energy range for all the nucleons in the state /; the cancelation between the 
3p2h and 2plh contributions of each Goldstone diagram does not occur at a 
given value of the energy of one of the nucleons contained in /. Note also that 
the numbers Nf^ and Nf^N') prevent the above cancelation to occur for 
the total numbers in Eqs. (pQ) too. Thus, QIT diagrams surely play a role in 
nucleon spectra. 

Let us emphasize some consequences of the sum rule condition in Eq. ([3]). 
First, while T n and T p are not observable, the total decay width Tnm (which is 
unaffected by FSI) is an observable. Second, it is obvious that the cancellation 
in Eqs. (T2]) and occurs because some terms are negative: these terms are 
necessarily QIT. It is well known that FSI are important; in addition the 
second property above tells us that the magnitude of QIT is the same as 
the one of plain FSI in the evaluation of T^m- This proportion can change 
in the evaluation of the spectra due to the presence of the numbers iV/(;v) 
and Nf(NN'), but in no way QIT can be neglected: without QIT, the emission 
spectra would be overestimated. It is important to stress that this discussion 
refers to a quantum mechanical treatment only. 

Before discussing our numerical results, it should be mentioned that the val- 
ues of f jv, f NN i and F it f are obtained from the corresponding A self-energy 
Goldstone diagrams. Explicit expressions for diagrams D and E in Fig. [T] are 
found in [13]. For the formal derivation of the decay rates from the pp, ph, 
hh, pp', ph' and hh' Feynman diagrams in the same figure, we proceed as fol- 
lows. The evaluation of the pp, ph and hh Goldstone diagrams contributing 
to 2N decays can be found in [TJ] (pp in the main text, the other two in the 
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Appendix). The expressions for the same time-orderings for pp' , ph! and hh' 
diagrams are shown in the Appendix of [15]. The remaining time-ordering for- 
mulas for pp, ph, hh, pp', ph' and hh' can be obtained from the just mentioned 
ones after some work: in any case, the main ingredients are already present in 
the mentioned works. 

As a starting point for the numerical analysis, in Table [1] we show the decay 
rates obtained for \ 2 C. We also predict that r np /T 2 = 0.84, r pp /T 2 = 0.12 
and r nn /r 2 = 0.04. The agreement with the recent KEK and FINUDA data 
is quite satisfactory. 

Table 1 

The non-mesonic decay widths predicted for ^ 2 C (in units of the free decay rate). 
The most recent data, from KEK-E508 [16J and FINUDA [T7], are also given. 

Decay rate Our KEK-E508 FINUDA 



T n 0.15 0.23 ±0.08 

T p 0.47 0.45 ±0.10 

Ti 0.62 0.68 ±0.13 

T 2 0.36 0.27 ±0.13 

T NM 0.99 0.95 ± 0.04 

r n /T p 0.33 0.51 ±0.13 ±0.05 

r 2 /r N M 0.37 0.29 ±0.13 0.24 ±0.10 



We turn now to our main concern in this Letter, which is the study of the 
emission spectra. The obtained emission spectra depend mainly on the phase 
space and the nucleon-nucleon strong potential. The fact that these spectra are 
normalized per non-mesonic weak decay makes the dependence on the weak 
interaction to be very small. In the present contribution we have employed 
the weak and strong interaction models used in [5]. These parameterizations 
turned out to reproduce the full set of decay widths for carbon hypernuclei 
determined recently at KEK |16j . Note that the use of the same strong inter- 
action model for both GSC and FSI not only provides internal consistency to 
the model, but gives us a criterion to select the strong interaction potential, 
since modifications of this potential leads to testable changes in both r 2 /r NM 
and the nucleon spectra. 

In Fig. [3] we show the neutron and proton kinetic energy spectra for the non- 
mesonic decay of \ 2 C. The dashed curves are the distributions of the IN decay 
nucleons (normalized per IN decay): as expected, they show a maximum at 
half of the Q-value for ^ 2 C non-mesonic decay and a bell-type shape due to the 
nucleon Fermi motion and the A momentum distribution in the hypernucleus. 
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The inclusion of 2N and FSI-induced decay processes provides the results 
given by the continuous lines (normalized per non-mesonic decay) and leads 
to a reduction of the nucleon average energy, thus filling the low-energy part of 
the spectrum and emptying the high-energy region. This outcome is explained 
as follows. The dominant contribution to FSI is provided by diagram Ufa): 
according to phase space arguments, when a cut on 3p2h states is considered 
one obtains a distribution that decreases monotonically with T/v, while a cut 
on 2plh states provides a bell-shaped, negative QIT peaked at T N ~ 70 MeV. 
Another important contribution is that of diagram Efb), which admits ?>p2h 
final states only: it produces positive and monotonically decreasing nucleon 
distributions with the same order of magnitude of the ?>p2h term of diagram 
Efa). While we reproduce fairly well the KEK neutron spectra, a rather strong 
overestimation is found of the KEK-E369 and FINUDA proton distributions. 
Our proton spectrum is instead closer to the old BNL-KEK data. 




T n (MeV) T p (MeV) 

Fig. 3. Neutron and proton kinetic energy spectra for j^ 2 C non-mesonic weak de- 
cay. The dashed (continuous, dot-dashed) lines are normalized per IN decay (per 
non-mesonic decay). Experimental data are from KEK-E369 [18], KEK-E508 [19] . 
FINUDA P2] and BNL-CERN [20]. 



Also from Fig. [3] we see a certain dispersion among the various experimental 
proton spectra (for instance, FINUDA shows a peaking structure at T p ~ 80 
MeV which is not seen in the other data nor in this and in previous calcu- 
lations). Compared with the nucleon spectra obtained with the INC model 
for FSI in [6] (dot-dashed curves), we obtain here similar results for T/v > 40 
MeV; the distributions of [6] for smaller T/v increase well beyond the predic- 
tions of the present work, which are more reliable than the ones based on the 
INC in this energy regime. Given the obvious differences between the present 
microscopic approach and the INC rescattering model, the previous agreement 
is a remarkable result. The QIT of the present approach replace the nucleon 
rescattering present in the INC, whose effect is to reduce the energy of the 
nucleons propagating through the residual nucleus. 

The opening angle distributions of nn and np pairs are reported in Fig. [4] To 
adhere to the KEK data, the predictions of the full model are obtained for a 
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30 MeV nucleoli kinetic energy threshold T^ 1 . The distributions from the IN 
decay (dashed curves) are strongly peaked at 6nn' = 180°. The QIT again 
have a crucial effect: they considerably reduce the back-to-back contribution, 
thanks mainly to diagram [2]^a) with 2plh final states; the non back-to-back 
region is strongly populated mainly by the 3p2/i contributions of diagrams [2]^a) 
and [2](b) (see continuous curves). The final results for the angular correlations 
turn out to be somewhat less back-to-back peaked than what found in [6] 
(dot-dashed curves), and, as in that work, are very sensitive to the value 
adopted for Tjy . The agreement with KEK-E508 data is rather good for the nn 
spectrum, while for np pairs a significant overestimation is obtained. Note that 
we reproduce for the first time the experimental nn opening angle spectrum: 
it is the quantum-mechanical nature of our scheme, i.e., the relevance of QIT, 
which brings to this achievement. Our overestimation of the KEK-E508 np 
distribution is compatible with the overestimation of the proton spectrum 
obtained in the same experiment (see Fig. [3]). 
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Fig. 4. Opening angle distribution of nn and np pairs. Normalization is as in Fig. [3j 
Data are from KEK-E508 [21]. 



In Fig. [5] we give the two-nucleon momentum correlation spectra, i.e., the nn 
and np distributions as a function of the momentum sum Pnn 1 = \pn +Pn'\ of 
two of the outgoing nucleons. The dashed lines correspond to the IN decay, 
while the continuous curves show the full result, with IN, 2N and FSI included 
together with a nucleon kinetic energy threshold T^ h = 30 MeV, as in the 
data also shown in the figures. As noted in Ref . [16] , the minimum in both the 
nn and np KEK-E508 distributions is mainly due to the low statistics and 
detection efficiency for events with pnn' > 350 MeV/c (the KEK detector 
geometry being optimized for back-to-back coincidence events); indeed, such 
dip structure has not been found in our calculation, which overestimates the 
data for large correlation momenta (especially for N np , consistently with the 
spectra discussed so far). 

The distributions of Fig. [5] at low momentum sum (say below 400 MeV/c) 
are mainly due to IN decays (which are strongly back-to-back correlated), 
while for higher momenta the contribution of 2N and FSI-induced decays is 
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p nn (MeV/c) P np (MeV/c) 

Fig. 5. Momentum correlation spectra of nn and np pairs, with pnn' = \pn ~\~Pn'\- 
Normalization is as in Fig. [3l Data are from KEK-E508 [16J . 

dominant (and produces less back-to-back correlated pairs). This behavior is 
confirmed by the momentum correlation of the sum N nn + N np shown in Fig. [6] 
for the opening angle regions with cos^/va" < —0.7 and cosOnn' > —0.7. 
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Fig. 6. Momentum correlation spectra for the sum of the nn and np pair numbers 
for the back-to-back (cosOnn < —0.7, continuous line) and non back-to-back 
kinematics (cosOnn > —0.7, dashed line). Normalization is per non-mesonic weak 
decay. Data are from KEK-E508 [22]. 



In Table [2] we show results for the coincidence numbers N nn , N np and N pp 
obtained for Tjy = 30 MeV and two angular regions, cos9 NN / < —0.7 and 
cos 8 A? n> > —0.7. Previous predictions of finite nucleus [H] and nuclear matter 
approaches [6], both based on the INC, and experimental data are also given. 
Again, a remarkable result is the fact that the present predictions are very 
similar to the ones of [6], while one notes less agreement with [IT], mainly 
due to the different models adopted to describe 2N decays. Comparison with 
data shows an overestimation of N np and N pp but rather good results for N nn . 
This confirms a systematic overestimation of the proton emission reported by 
KEK-E508. 



Summarizing, a microscopic approach including GSC and FSI on the same 
footing is used to evaluate the nucleon emission spectra in non-mesonic weak 
decay of hypernuclei. Within our microscopic model, QIT play a key role: in 
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Table 2 

The nucleon coincidence numbers are given for T^ h = 30 MeV and the angular 
regions with cos 9nn* < —0.7 and cos Onn 1 > —0.7 (in parenthesis). Data are from 
KEK-E508 [ZQ. 





N nn 


N np 


N pp 


This work 


0.11 (0.13) 


0.35 (0.31) 


0.04 (0.05) 


INC p] 


0.15 (0.18) 


0.35 (0.52) 


0.08 (0.27) 


INC [6] 


0.11 (0.10) 


0.30 (0.25) 


0.05 (0.07) 


KEK-E508 


0.083 ±0.014 


0.138 ±0.014 


0.005 ± 0.002 




(0.083 ± 0.020) 


(0.060 ± 0.018) 





the single-nucleon emission spectra they are responsible for moving intensity 
from the high-energy region to the low-energy region, while in the opening 
angle distribution of nucleon pairs the strong reduction of the back-to-back 
peak is entirely due to QIT. Discrepancies with experiment remain, but are 
relegated to spectra involving protons. Further work is in order to understand 
such a disagreement. A forthcoming coincidence experiment at J-PARC [23] 
will allow a measurement of the nucleon spectra and a determination of T n , T p 
and T2 for \ 2 C with improved accuracy. Moreover, our microscopic approach is 
particularly suitable for the inclusion of new decay mechanisms introduced by 
fermion antisymmetrization and the A(1232)-resonance. These contributions 
deserve consideration. 
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